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M.  ABSTRACT 
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A detailed,  human  ballistic  mortality  model  has  been  developed  at  the 
Ballistic  Research  Laboratory  which  considers  the  interdependence  of  the  dif- 
ferent organs  and  tissues  in  a human.  Variations  by  individuals  in  their 
vulnerabilities  to  ballistic  damage  are  described  by  probability  density 
functions.  A Monte  Carlo  solution  of  the  model  for  calculating  the  mortality 
rate  of  a large  number  of  soldiers  exposed  to  a common  ballistic  threat  is 
outlined.  Representative  wounds  are  constructed  in  a human  Phantom  in  v_^ 
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I.  INTRODUCTION 


A methodology  which  can  be  used  to  calculate  the  incapacitation  rate 
of  soldiers  exposed  to  a ballistic  threat  is  needed  by  the  Army.  These 
rates  would  be  used  by  military  planners  to  develop  more  efficient  offen- 
sive weapons  and  to  optimize  the  development  of  defensive  material  such 
as  body  armor.  The  results  could  also  be  used  in  the  planning  of  both 
offensive  and  defensive  tactics. 


A statistical  model  which  predicts  the  mortality  rate  of  a large 
number  of  soldiers  exposed  to  a ballistic  threat  has  been  developed  by 
Beverly.  This  model  assumes  that  the  mortality  rate  of  a large  number 
of  soldiers  in  battle  can  be  modeled  by  the  phantom  description  of  an 
average  soldier.  The  Phantom  will  serve  as  the  target  for  ballistic 
threats  in  computerized  enactments  of  battle  scenes.  Representative 
projectiles  are  picked  from  the  threat  and  projected  toward  the  Phantom. 
A representative  wound  is  constructed  for  each  projectile  that  striker 
the  Phantom.  The  mortality  probability  is  calculated  for  each  wound  and 
the  mean  probability  is  calculated  for  a large  number  of  wounds  to  oo- 
tain  the  mortality  rate  expected  from  the  threat. 


The  vulnei  ibility  (VI)  to  ballistic  injury*  of  an  organ  in  an 
individual  sold.sr  is  defined  as  the  average  threshold  volume  of  tissue 
destruction  in  t ie  organ  that  would  lead  to  his  death.  This  definition 
implies  that  the  soldier  has  suffered  no  injury  to  his  other  organs. 

The  vulnerability  of  the  same  organ  in  the  Phantom  is  described  by  a 
probability  density  function  (PDF)  which  predicts  the  frequency  of  oc- 
currence of  VI -values  in  the  total  soldier  population. 


by 


A damage  Di  is  defined  for  an  injury  to  the  ith  organ  in  a soldier 


D. 

1 


(1) 


Where 

WL  = the  volume  of  tissue  destruction  i.n  the  ith  organ. 

VI.  = the  mortal  threshold  volume  of  destruction  in  his 
ith  organ 


An  addition  law  is  assumed  to  exist  which  can  be  used  to  accumulate 
the  damages  suffered  by  different  organs  to  obtain  a damage  D for  the 
total  wound.  This  addition  law  is  averaged  over  the  total  soldier 


1. 


W.B.  Beverly,  " A Human  Ballistic  Mortality  Model,"  to  be  published  as 
as  a BEL  report. 


population  to  obtain  the  law  which  will  be  used  in  Phantom  mortality  cal- 
culations. This  rule  is  assumed  to  be  expressed  analytically  as 


D - FCD^ 


(2A) 


A wound  in  the  Phantom  whose  damage  is  calculated  to  be  1 or  greater 
is  assumed  to  be  mortal. 


Baker  et  al2'3  have  conducted  studies  which  indicate  the  function 
F(.D^)  can  be  approximated  by  the  quadratic  addition  of  organ  damages. 

Equation  2A  could  then  be  expressed  as 


[*  °i2] 


1/2 


(2B) 


where 


I » the  number  of  organs  in  a soldier 


i * the  running  index  over  organs 

We  will  use  the  quadratic  rule  for  the  present  but  further  studies 
should  be  conducted  to  verify  its  accuracy. 


The  mortality  probability  of  the  Phantom  (or  the  mortality  rate  of 
the  soldiers)  for  a particular  wound  is  the  product  of  the  probability 
of  the  wound  occurring  and  the  death  rate  which  would  result  from  the 
wound.  The  mortality  probability  M of  the  Phantom,  when  he  is  exposed 
to  a ballistic  threat,  is  the  average  of  his  mortality  probabilities  for 
all  possible  wounds  and  is  given  in  reference  1 as 


M -f  j BT(S(rj,p)-^] j H(D-l)$wfD)  dDdV  dP  (3A) 

L'J  *l(Vil),,-*I(VII)/6lD"F(Di)J<1V2  d(VI1)...d(VlI)  (3B) 


’ S.P . Baker,  B.  O’Neil,  W,  Hadden,  Jr.,  and  W.B.  Long,  "The  Injury 
Severity  Score:  A Method  for  Describing  Patients  with  Multiple 
Injuries  and  Evaluating  Emergency  Care , " THE  JOURNAL  OF  TRAUMA,  VOL. 
14,  No.  3,  1974. 

S.P . Baker  and  B.  O’Neil,  "The  Injury  Severity  Score:  An  Update ,” 

THE  JOURNAL  OF  TRAUMA,  VOL  16,  No.  11,  1976. 
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where 


rj  ■ the  origin  of  the  projectile  in  the  ballistic  threat  coordinate 
system 

dVj  « the  infinitesimal  volume  in  the  ballistic  threat  coordinate 
system 

p » the  momentum  of  the  projectile 

dP  * the  infinitesimal  volume  in  projectile  momentum  space 

w * the  wound  volume  in  the  Phantom 

s(*l»P)  a the  projectile  source  term  normalized  to  one  source 
projectile 

BT[S(rlfp)->w3  ■ the  function  that  describes  the  production  of  a 
wound  by  a projectile 

H(D-l)  * the  Heaviside  step  function 

*i(VI.)  * the  probability  density  function  which  predicts  the  mortal 
J volumes  of  tissue  destruction  VI.  for  the  i**1  organ  of 
soldiers  1 

5 [D-F(D^)]  = the  Kronecker  delta  function 

dV2  = the  infinitesimal  volume  in  the  Phantom  coordinate  system. 

It  should  be  noted  that  the  innermost  integral  of  equation  3B  is  used  to 
calculate  the  organ  destruction  volumes  VWj_  which  are  used  to  calculate 
the  Dj.  The  quantity  D is  constant  insofar  as  this  integration  is 
concerned  and  is  used  in  the  Kronecker  delta  functional  to  evaluate  $ 
as  a function  of  D. 

An  analytic  solution  to  equation  3 would  usually  require  many 
cumbersome  calculations.  Such  a heroic  effort  can  be  avoided  by  develop- 
ing a completely  stochastic  solution  where  all  calculations  are  relatively 
simple.  However,  the  calculations  will  need  to  be  repeated  many  times 
until  estimates  have  been  calculated  to  the  desired  statistical  accuracy. 
The  proposed  Monte  Carlo  solution  to  the  mortality  model  is  briefly 
outlined  here  and  will  be  discussed  in  detail  in  the  next  section  of 
this  report. 

Representative  projectiles  are  picked  from  the  PDF's  that  describe 
the  ballistic  threat  and  are  projected  toward  the  Phantom.  Representative 
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wounds,  using  the  Dubin  Cavity  Model,4  are  constructed  for  those  projec- 
tiles which  strike  the  Phantom.  An  estimate  of  the  volume  of  destruction 
inflicted  upon  each  organ  by  a wound  is  obtained  from  the  number  of 
sample  prints,  described  later,  lying  within  that  organ.  A threshold 
mortality  volume  VI^  is  picked  for  each  injured  organ  by  sampling  the 

appropriate  PDF  and  is  used  to  calculate  the  damage  inflicted  upon  that 
organ.  The  total  damage  of  the  wound  is  calculated  using  equations 
2A  and  2B.  A kill  scored  for  a wound  when  its  damage  is  one  or  greater 
and  a survival  is  scored  when  its  damage  is  less  than  one.  'Hie  effects 
of  a large  number  of  similarly  derived  wounds  are  evaluated  and  averaged 
and  the  mortality  rate  of  soldiers  is  calculated  as  the  expected  mortal- 
ity probability  of  the  Phantom. 


IX.  THE  STOCHASTIC  SOLUTION 

Section  II  of  this  report  is  divided  into  three  sub-sections  where 
the  first  sub-section  is  used  to  describ*  the  general  procedures  used  to 
calculate  the  mortality  score  of  a wouna.  These  procedures  are  discussed 
as  if  the  step-by-step  operation  of  a hypothetical  computer  program  were 
being  described.  The  Monte  Carlo  sampling  techniques  that  are  used  here 
to  pick  VI -values  for  injured  organs  are  described  in  detail  in  sub- 
section IB.  The  Monte  Carlo  techniques  that  are  used  here  to  estimate  the 
volumes  of  destroyed  tissue  in  the  injured  organs  are  described  in  detail 
in  sub-section  C. 

A.  The  General  Solution 

■> 

A projectile  having  momentum  p ...s  picked  at  origin  r. , by  sampling 
the  PDF  which  describes  the  ballistic  threat.  The  picking  and  tracking 
of  projectiles  which  do  not  intercept  the  Phantom  is  reduced  by  picking 
only  projectiles  whose  momentum  vectors  lie  within  the  solid  angle  sub- 
tended by  the  Phantom- enclosing- sphere  (see  Figure  1) , 5 The  projectile- 
parameter  picking  teheniques  are  not  discussed  in  this  study  except  to 
assume  that  the  threat  description  can  be  efficiently  sampled  using 
standard  sampling  techniques. 

A shotline  is  projected  into  the  Phantom  for  each  projectile  that 
intercepts  the  Phantom.  The  deflection,  break-up,  and/or  deformation  of 
projectiles  which  strike  bone  tissue  are  not  included  in  this  study  but 
should  be  ultimately  considered.  The  present  discussion  will  proceed  as 
if  the  shctline  of  each  projectile  were  a straight  line. 

H. r.  Dubin , "A  Cavitation  Model  for  Kinetio  Energy  Projectiles  Pene- 
trating Gelatin,  " Ballistic  Research  Laboratories  Memorandum  Report 
So.  2422,  December  1974.  (AD  #B000910L) 

5'W.B.  Beverly,  "RAYMAN:  A FORTRAN  Computer  Code  for  Tracing  Rays 
Through  a Detailed  Human  Phantom,"  Ballistic  Research  Laboratory 
Technical  Report  ARBRL-TR-2030,  November  1977 . (AD  #A051057) 


The  point  of  entry  and  exit  of  a shot  line  into  different  organs 
along  its  path  are  located.  The  existing  GIFT* computer  program  can  be 
used  to  perforin  this  calculation  when  the  Phantom  is  described  using  the 
combinatorial  geometry  (COM-GEOM)  technique.  The  geometry  and  input  re- 
quirements used  in  GIFT  are  discussed  in  detail  in  the  G.tFT  User's  Manual.7 

A shotline  within  an  organ  is  divided  into  segments  where  all  seg- 
ments, except  possibly  the  last  one  which  includes  the  exit  point,  have 
equal  lengths.  A truncated  right  cone  (TRC)  is  constructed  about  each 
segment  to  represent  the  wound.  The  radii  of  the  cones  are  calculated 
using  techniquea  developed  in  Reference  4.  The  radii  of  a cross  section 
of  a TRC  will,  in  general,  decrease  with  increasing  penetration  within 
an  organ,  Each  cone  is  constructed  using  the  physical  characteristics 
of  the  tissue  type  (primary  tissue)  that  is  located  on  its  axis.  The 
cross  section  of  a typical  wound  construction  is  displayed  in  Figure  2. 

Points  are  picked  with  equal  probability  over  the  volume  of  a wound 
by  conducting  a sampling  pass,  in  sequence,  over  each  TRC  within  the 
wound.  We  assume  that,  on  the  average,  one  point  is  picked  from  each 
unit  of  volume.  A score  of  1 is  assigned  to  each  point  which  is  em- 
bedded within  the  primary  organ  and  this  score  is  added  to  a bin  reserved 
for  that  organ. 

A point  which  lies  outside  the  primary  organ  requires  additional 
treatment  as  the  energy  flux  reaching  the  point  may  not  be  sufficient  to 
destroy  the  tissue.  A projectile  is  assumed  to  not  intercept  the 
secondary  organ  at  any  poi.it  within  the  TRC  being  sampled.  Some  of  the 
energy  projected  toward  the  secondary  organ  is  reflected  at  the  inter- 
face. This  reflection  results  in  a lower  energy- flux  level  in  the 
secondary  organ  than  would  be  the  case  if  no  interface  were  present. 
Additionally,  the  energy-flux  level  required  to  destroy  tissue  in  the 
secondary  organ  can  differ  from  that  required  in  the  primary  organ. 

These  effects  will  cause  the  wound  cavity  in  the  secondary  organ  tc 
differ  from  the  cavity  predicted  for  homogeneous  tissue.  The  effect  is 
shown  in  Figure  3 as  an  expanded  view  of  Region  B in  Figure  2. 

A radius  vector  is  constructed  from  the  sample  point  in  the  secon- 
dary organ  to  the  axis  of  the  TRC.  The  exit  point  from  the  primary 


ft 

‘ W.  Guber , R.  Nagel,  P.  Goldstein , P.S.  Mittleman,  and  M.  Kalos,  "A 
Geometric  Description  Technique  Suitable  for  Computer  Analysis  of  Both 
the  Nuclear  and  Conventional  Vulnerability  of  Armored  Military 
Vehicles,"  Mathematical  Applications  Group,  Inc.,  Report  No.  MAGI  6701, 
August  1967. 

7 

'L.A.  Bains,  Jr.  and  M.J.  Reisinger,  " The  GIFT  Code  User  Manual ; VOL  I, 
Introduction  and  Input  Requirements,  " Ballistic  Research  Laboratory 
Report  No.  1802,  July  1975 . (AD  #B006037L) 
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organ  along  this  ray  is  than  located.  Ilia  ana rgy-ra fleet ion  coefficient 
of  the  two  organs  is  used  to  calculate  the  energy  flux  transacted  across 
the  interface.  The  energy  flux  transacted  to  the  saaple  point  is  cal- 
culated and  a score  of  1 is  tallied  for  the  secondary  organ  when  the  flux 
is  sufficient  to  destory  the  secondary  tissue.  A zero  score  is  tallied 
when  the  energy  flux  arriving  at  the  point  is  not  sufficient  to  destroy 
the  tissue. 

A stochastic  estimate  of  the  volute  of  tissue  destroyed  in  each 
organ  is  calculated  using  the  results  of  several  sampling  passes  over 
a wound  voluae.  Assuaing  that  one  point  was  picked  per  pass  from  each 
unit  of  volume,  an  estinate  of  the  volimw  of  tissue  destroyed  in  the  ith 
organ  may  be  calculated  using 

wi  * r£vik  (4) 

where 

k - the  running  index  over  sampling  passes 
K = the  number  of  sailing  passes 

V 

ik  * the  total  nuaber  of  points  (scores)  accumulated 
for  the  ith  organ  during  the  kth  pass. 

The  standard  deviation  «(VW  ) of  this  estimate  is  calculated  using  8 


The  bar  is  placed  over  a variable  to  show  that  a mean  estimate  of  its 
value  is  being  used. 

A value  of  of  each  injured  organ  is  picked  by  sampling  its 

probability  density  function.  The  damage  to  each  injured  organ 
(equation  1)  and  the  total  damage  (equation  2B)  are  calculated  yielding 


D * 


(6A) 


8. 


I.  Beers,  "Introduction  to  the  Theory  of  Errors,  " Addison-Wesley 
Publishing  Co.,  Inc.,  1953. 
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The  variances  of  these  damage  estimates  are  derived  to  be8 


52D  * 


1 

D* 


(6B) 


A criterion  is  set  for  the  minimum  acceptable  statistical  uncertainty 
in  D.  Another  series  of  sampling  sweeps  will  be  conducted  when  this 
criterion  is  not  reached  during  the  current  series.  A kill  is  then 
scored  for  a wound  whose  daaage  estimate  is  one  or  greater. 

As  noted  earlier  in  this  section,  only  those  projectiles  which 
have  a high  probability  of  striking  the  Phantom  are  picked  and  tracked. 
The  mortality  rate  due  to  a threat  whose  projectiles  are  emitted  iso- 
tropically about  their  origin  points  is  given  by 

"■?J H<D„  - » (7> 

where 

n = the  running  index  over  projectile  histories 

N = the  number  of  projectile  histories 

Dn  * the  total  damage  of  the  n^  wound 

0 (rt)  = the  solid  angle  subtended  by  the  Phantom-enclosing 
sphere  at  the  projectile  origin 

H(Dn  - 1)  * Heaviside  step  function. 


The  normalization  cf  the  mortality  rate  for  an  anisotropic  source  term 
is  not  discussed  in  the  present  study. 

The  general  logic  used  by  the  hypothetical  computer  program  to  cal- 
culate an  estimate  of  the  mortality  rate  due  to  a ballistic  threat  is  now 
outlined.  Detailed  discussions  of  the  Monte  Carlo  techniques  used  here 
are  given  in  Sections  II. B and  II. C. 

1.  Pick  a projectile  having  momentum  p at  origin  by  sampling 
the  ballistic  threat  probability  density  functions. 

2.  Track  the  projectile  through  the  Phantom-enclosing-sphere 
and  construct  a representative  woiud  if  the  projectile  strikes  the 
Phantom. 


3.  Zero  the  score-hod  array  HOD,  the  score  accumulation  array 
SCORE,  and  the  squared-score  accumulation  array  SC0RE2.  The  array  HOD 
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is  used  to  accumulate  the  scores  acquired  by  each  organ  during  a pass. 
The  arrays  SCORE  and  SCORE 2 are  used  to  accumulate  the  scores  and 
squared  scores  gathered  over  several  passes. 
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4.  Conduct  a stapling  pass  over  the  wound  where  points  are  picked 
with  equal  probability  over  the  volume  of  the  wound.  Determine  the  organ 
in  which  each  point  is  embedded. 

5.  Add  the  score  of  each  sample  point  to  the  HOD  bin  reserved  for 
the  organ  in  which  the  point  is  embedded.  A non- zero  score  is  the 
reciprocal  of  sample-point  density.  This  definition  of  score  value 
yielded  the  scores  of  1 and  0 which  were  used  earlier  in  this  report.  A 
sample-point  density  of  one  point  per  unit  volume  was  used  in  this  case. 

6.  At  the  completion  of  each  pass,  add  the  scores  accumulated  in 
HOD  to  SCORE.  .Square  each  score  in  HOD  and  add  the  result  to  SC0RE2.  A 
bin  is  reserved  in  SCORE  and  SC0RE2  for  each  injured  organ. 


7.  Zero  all  bins  in  HOD  and  conduct  K-l  similar  score  passes. 


8.  Calculate  an  estimate  of  the  volume  of  tissue  destruction 
using  the  current  version  of  equation  4, 


SCORE (i) 

ic 


(8A) 


9.  Calculate  an  estimate  of  the  variance  62(VWj)  of  each  volume 
estimate  using  the  current  version  of  equation  5, 


SC0RE2 (i)  - K 
62(VW.)  " K(K  - X) 


VW 


(8B) 


10.  Pick  a mortal  volume  VI  for  each  injured  organ  by  sampling 

the  appropriate  PDF.  Calculate  an  estimate  of  the  damage  (equation  1) 
and  the  standard  deviation  of  the  damage  estimate  for  each  injured 
organ.  The  damage-estimate  variance  is  given  by 


6 2D  - 


2 


(SC) 


1 

1 


i 


1 


\ 

i 


l 

■1 

3 


i 

4 

\ 


11.  Calculate  an  estimate  of  the  total  damage  D of  the  wound 
(equation  2B)  and  the  standard  deviation  of  the  damage  estimate 
(equation  6B) . 

12.  Compare  the  variance  of  the  damage  estimate  (equation  8C) 
with  the  maximum  variance  threshold  criterion.  Proceed  to  the  next 
step  when  the  variance  is  less  than  this  threshold  value.  Conduct 
another  series  of  sweeps  over  the  wound  when  the  variance  is  too  large. 


j 

i 


i 

i 


y? 
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13.  A kill  is  scored  when  the  damage  equals  or  exceeds  the  mortal 
damage.  The  geometrical  normalization  is  calculated  for  a killing  pro- 
jectile using  equation  7. 

14.  Conduct  similar  calculations  until N projectiles  histories 
have  been  completed.  The  killing  probability  of  the  threat  is  then 
calculated  using  equation  7. 

B.  Sampling  the  VI  Probability  Density  Functions 

The  probability  density  function  of  the  mortal  volumes  of  tissue 
destruction  is  measured  or  calculated  as  a histogram.  These  functions 
can  be  derived  by  an  analysis  of  gunshot  case  histories  in  combination 
with  assessments  by  shock-trauma  doctors  of  wound  effects.  Assuming 
the  histogramic  form,  Sample  VI  values  can  be  picked  using  the  funda- 
mental principle  of  Monte  Carlo  (FPMC).9  The  application  of  this 
principle  yields  a sample  value  given  by 

fVIi 

I MVI)d(VI)  = RN(0,1)  (9) 

'o 

where 

RN(0,1)  = a random  number  picked  with  equal 
probability  on  the  range  of  0 to  1 

= the  picked  VI  value  for  the  ith  organ 
Each  4^  (VI)  is  normalized  to  one  over  its  complete  range  of  values. 

The  sampling  process  is  hastened  by  evaluating  the  integral  in 
equation  (9)  at  each  histogramic  bin  interval  and  storing  the  results  in 
a VI  versus  CDF  (cumulative  distribution  function)  table.  These  values 
are  assumed  to  be  stored  in  ascending  order  where  the  last  CDF -value 
in  each  table  is  one.  The  plot  of  a typical  PDF  and  its  associated  CDF 
are  displayed  in  Figure  4. 

A sample  VT-value  may  now  be  picked  using  a procedure  where  the  only 
operations  are  the  comparison  of  the  table  CDF-values  with  a random 
number  followed  by  a linear  interpolation.  This  procedure  is  illustrated 
in  Figure  4 and  is  outlined  below.  We  will  assume  that  the  VI  versus  CDF 
values  of  the  organ  in  this  example  are  stored  in  the  CUM  array.  It 
should  be  noted  that  J and  K which  are  used  here  apply  only  to  equation 
10  and  this  definition  should  not  be  used  elsewhere  in  this  report.  The 
VI  values  are  stored  in  the  top  row  of  the  CUM  matrix  (array)  and  the  CDF 
are  stored  in  the  bottom  row  of  the  CUM  array. 


9. 


y.A,  Sahreider,  " The  Mo. its  Carlo  Method »"  Perganion  Press,  1966. 
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1.  Calculate  a random  number  RN(0,1). 

2.  Lcoate  the  adjacent  pair  of  COP  values  in  CUM  which  bound 
RN(0, 1) . The  smallest  value  is  indexed  by  J and  the  largest  value  is 
indexed  by  K. 


3.  The  sample  Vl-value  is  calculated  using 


(10) 


C.  The  Estimation  of  Organ  Destruction  Volumes 

The  volume  of  tissue  destroyed  in  each  organ  of  the  Phantom  is  also 
estimated  stochastically.  This  estimate  is  calculated  by  picking  points 
with  equal  probability  over  the  volume  of  the  wound  and  identifying  the 
organs  in  which  they  are  embedded.  An  estimate  of  the  volume  of  tissue 
destroyed  in  a particular  organ  is  proportional  to  the  number  of  non- 
zero score  points  that  are  located  in  the  organ.  The  statistical  ac- 
curacy of  this  estimate  depends  upon  the  density  of  sampling  points 
oicked  within  the  wound  and  will  converge  toward  the  true  value  as  the 
density  of  sample  points  becomes  large. 

A wound  is  initially  constructed  as  a series  of  coaxial,  truncated 
right  cones.  The  conical  geometry  is  modified  at  m organ  interface  by 
the  albedo  correction  but  sampling  points  are  still  picked  with  equal 
probability  over  the  volume  of  the  original  cone.  However,  a zero  score 
is  tallied  for  points  that  lie  outside  the  modified  wound  volume. 

A TRC  is  divided  into  two  parts  to  facilitate  the  calculation  of 
sample  points  (see  Figure  5).  A coaxial  cylinder  (RCC),  which  has  a 
radius  equal  to  the  radius  at  the  small  end  of  the  TRC,  is  constructed 
in  the  TRC.  Sampling  points  can  be  easily  picked  over  this  volume  by 
using  the  FPMC.  The  picking  of  sample  points  in  the  remainder  of  a cone 
(REM)  will  require  the  use  of  the  less  efficient  rejection  technique. 

A sampling  pass  over  a TRC  will  be  first  conducted  over  its  RCC  and  then 
over  its  REM.  The  calculation  of  an  organ  destruction  volume  estimate 
is  now  discussed  in  step-by-step  detail. 

1.  Calculate  the  volume  of  the  TRC  and  identify  this  quantity  as 
VCN.  Calculate  the  volume  of  the  enclosed  cylinder  and  identify  this 
quantity  as  VCY.  Truncate  VCN  to  an  integar  and  identify  this  number  as 
M.  This  definition  of  M is  only  used  in  this  sequence  of  steps  and 
should  not  be  confused  with  the  mortality  rate  M used  earlier  in  this 
report.  Start  an  associated  running  index  m at  1. 


2.  Calculate  a random  number  RN^(0,1).  The  subscript  is  added 
to  the  usual  identification  of  a random  number  to  identify  its  origin  for 
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those  cases  where  a particular  random  number  may  be  translated  to  a new 
range  and  reused.  All  unsubscripted  random  numbers  are  used  once  and 
discarded. 

3.  Compare  RNj(m-l,»D)  with  VCY.  Go  to  step  1 of  the  F.EM  sampling 

if  RN(m-l,rc)  is  greater  than  VCY.  An  outline  of  REM  sampling  immediately 
follows  the  outline  of  RCC  sampling.  We  have  assumed  that  the  sampling 
density  is  one  point  per  unit  volume. 


4.  The  PDF  of  points  in  a volume  as  a function  cf  one  of  the 
coordinates  is  proportional  to  the  cross-sectional  area  of  the  volume 
as  a function  of  the  same  coordinate.  This  concept,  when  applied  to 
the  RCC  enclosed  by  a TRC,  permits  the  FPMC  to  be  used  to  calculate  the 
penetration  hm  of  a quadruplet  of  sample  points  by  solving 


RN(m-l.m) 


(11) 


where 


h * the  penetration  along  the  axis  of  the  RCC  as 
measured  from  the  small  end  of  the  associated 
TRC 

R * the  radius  of  the  cylinder. 

S.  The  PDF  of  points  on  a plane  area  as  a function  of  one  cf 
the  coordinates  is  proportional  to  the  length  of  the  curve  within  the 
area  defined  by  the  same  coordinate.  This  concept  permits  the  FPMC 
to  be  used  to  calculate  the  radius  rffi  of  the  quadruplet  of  sample  points, 

located  on  the  cross  section  in  the  RCC  which  is  defined  by  h , by 
solving  m 


/ ^Trrd  = irR2  • RN(0,1)  (12) 

Jo 

where 


r * the  radius  of  a concentric  circle  within  the 
specified  circular  cross  section  of  the  RCC. 


6.  Locate  an  azimuthal  reference  on  the  circle  that  has  been 
defined  in  steps  4 and  5.  A sample  point  is  picked  with  equal  proba- 
bility on  each  90°  arc  of  the  circle  by  solving 


TT 

• 

2 


RN2(j-l,j) 


(13) 
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where 


d ■ the  azimuthal  angle  (radians) 
j ■ the  circle  quadrant  index. 

7.  Transfora  each  saaple  point  PCh^*  ? |#,  6^)  into  the  Phantom 

Cartesian  coordinate  system  to  obtain  four  points  (Pjm(x,y,z).  Determine 

the  organ  in  which  each  point  is  embedded.  The  sample  point  density  is 
four  points  per  unit  volume  so  the  non-zero  scores  are  now  0.25.  Add  a 
score  of  0.25  to  the  appropriate  HOD  bin  for  each  point  which  lies  within 
the  primary  organ.  Go  to  step  9 if  all  four  points  lie  within  the  pri- 
mary organ. 


8.  Calculate  the  energy  flux  which  is  transmitted  from  the  shot- 
line in  the  primary  point  in  the  secondary  organ.  Add  a score  of  0.25 

to  the  appropriate  bin  of  the  HOD  array  when  the  energy  flux  at  the  point 
exceeds  the  destruction  threshold  value.  A zero  score  is  added  to  the 
HOD  array  when  the  energy  flux  at  the  point  is  less  than  the  destruction 
threshold  value. 

9.  This  completes  the  sailing  in  the  first  unit  volume  of  the 
RCC.  Translate  RNjCm-l.m)  to  a new  range  using 

RNjCm.m+l)  *»  RN^(m-l,m)  * 1 . (14) 

Add  1 to  the  running  index  m and  return  to  step  3.  It  should  be  noted 
that  the  index  m is  incremented  after  the  translation  of  the  random 
number  so  that  RNi(m,m+l)  as  used  in  equation  14  is  equal  to  RN  ^m.m-l.m) 
as  used  in  step  3 when  the  return  is  completed. 

The  sampling  for  one  pass  through  the  RCC  is  completed  when  a jump 
to  the  first  step  of  REM-sampling  is  directed.  A sample  point  in  a REM 
is  also  chosen  by  picking,  in  order,  its  penetration  hffl,  radium  r^,  and 

azimuth  angles  9.  . The  coordinates  r uni  0.  are  calculated  using 
jm  m jm 

th^  same  techniques  that  weie  used  to  calculate  these  picks  during  RCC 
sampling.  However,  using  the  FPMC  to  pick  a penetration  hffi  would  re- 
quire finding  the  roots  of  a cubic  equation.  This  difficult  calculation 
is  avoided  by  using  the  rejection  technique  of  Reference  9 to  pick  the 
penetration  coordinate  of  sampled  points.  The  procedure  used  to  pick 
sample  points  in  a REM  is  outlined  below.  . 

1.  Compare  the  running  index  m with  M and  go  to  step  5 if  they 
are  equal. 


2.  The  PDF  of  points  in  a FEM  as  a function  of  the  penetration 
along  its  axis  (Figure  5)  is  given  by 
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where 


PDF(h)  * CA(h) 


(15A) 


h * the  penetration  along  the  axis  of  the  REM 
as  measured  from  the  small  end  of  the  TRC 

PDF (h)  * the  PDF  of  points  in  a REM 

C * the  normalization  constant. 


The  area  of  a cross  section  in  a REM  is  the  difference  between  the 
area  of  the  associated  cross  section  of  the  TRC  and  the  area  of  the 
enclosed  cross  section  of  the  RCC  and  is  given  by 


where 


A(h)  « irr2  - irR2  » Sh(2R  + Sh) 


(15B) 


S 


R2  - R 
- -fl 


(ISC) 


R « the  radius  of  the  TRC  at  its  small  end 
R2  - the  radius  of  the  TRC  at  its  large  end 
H « the  height  of  the  REB  (and  TRC) . 


The  procedure  for  picking  the  hm- coordinate  of  a sample  point  is 

illustrated  graphically  in  Figure  6 and  is  discussed  from  that  per 
spective  in  the  following  step-by-step  outline. 


2A.  Plot  the  function  PDF(h)  over  its  range  of  0 to  H. 

2B.  Construct  a rectangle  that  barely  enclosed  this  curve. 

2C.  Pick  a point  P(hm>PDFm)  at  anX  location  within  the 

rectangle  with  equal  probability.  The  coordinates  of  the  point  are 
calculated  using 


hm  = H • RN(0,1) 


(16A) 


PDFm  - PDF(H)  • RN(0,1)  (16B) 

2D.  The  quantity  hm  is  a good  choice  when  p(hm>pDpm)  lies  under 
the  curve  (PDF(h).  Proceed  to  the  next  step  when  this  favorable  event 
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Valid  Sample 
Invalid  Sample 


occurs.  Discard  the  current  point  P(hm,PDFJ  when  it  lies  outside  the 

n m 

curve  and  pick  a new  point  by  returning  to  step  2C.  Repeat  step  2C 
until  a fa.orable  point  is  picked. 

3.  Calculate  the  radius  coordinate  y_  of  the  sample  points 

using  B 

2m*dr  « ir£(R2)2  - R^2]  • RN(0,1)  (17A) 

R - R + h°i(Rg  “■&  (17B) 

m n 

The  reader  should  note  that  the  integration  interval  is  changed  in 
equation  (17A)  from  that  used  in  equation  (12)  for  RCC  sampling. 

4.  The  azimuthal  angle  coordinates  of  four  sample  points  are 
picked  using  the  same  techniques  that  were  outlined  in  the  RCC  sampling 
discussion.  Similarly,  the  embedding  organs  of  sample  points  are 
identified  and  scores  are  calculated  and  tallied.  The  running  index  is 
incremented  by  one  at  the  collet  ion  of  each  tallying  of  the  scores  for 
four  points  and  the  calculation  flow  is  returned  to  step  1 of  the  REM 
sampling. 

5.  A fractional  weight  W is  assigned  to  the  last  four  sample 
points  and  given  by 


W ■ VOI  - M (18) 

where  VON  is  the  TRC  volume  and  M is  the  largest  integer  in  VCN. 

The  coordinates  of  these  last  points  are  picked  using  steps  2-4  but 
the  scores  now  tallied  are  either  W/4  or  0. 

This  completes  the  outline  of  one  sampling  pass  through  a TRC. 

Each  TRC  in  a wound  will  be  similarly  sampled  during  a pass  and  an 
adequate  number  of  passes  will  be  conducted  over  the  wound  to  obtain 
damage  estimates  having  the  desired  statistical  accuracy. 


III.  PROGRESS  AND  PROGNOSIS 

A Monte  Carlo  solution  to  the  human  ballistic  mortality  model  has 
been  outlined.  Sampling  procedures  which  can  be  used  to  pick  variant 
values  are  developed  and  discussed.  Procedures  are  developed  for  cal- 
culating the  confidence  limits  on  damage  estimates.  Adequate  sampling 
of  a wound  volume  is  assured  by  comparing  the  variance  of  a damage  esti- 
mate with  the  maximum  acceptable  variance  criterion.  Further  sampling 
of  a wound  vclume  would  then  be  automatically  directed  by  the  computer 
program  until  the  variance  criterion  was  met. 
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The  probability  density  functions  that  predict  the  incapacitating 
volumes  of  destruction  for  the  different  organs  have  not  been  evaluated 
due  to  lack  of  experimental  data.  A data  bank  of  hospital  gunshot  cases 
which  could  be  used  for  this  purpose  is  being  collected  by  Sacco  and 
Sturdivan  but  the  acquisition  is  not  complete  at  this  time.  The 
precision  of  such  data  may  not  be  very  high  since  attending  physicians 
focus  their  attention  upon  saving  the  life  of  the  patients  rather  than 
making  accurate  experimental  measurements.  Mortality  assessments  of 
wounds  by  shock  trauma  doctors  have  been  collected  by  Kokir.akis  and 
Cooper lf but  the  latter  data  is  not  adequate  within  itself  to  derive 

the  necessary- yroeahi Ht.y  density  functions.  The  best  set  jaf  proba- 

bility  density  functions  will  probably  be  obtained  by  utilizing  both 
data  sources. 

Many  of  the  subroutines  and  much  of  the  logic  required  by  a 
mortality  computer  code  exist  in  other  codes.1*  »5*®  These  subroutines 
and  logic  will  need  to  be  assimilated  into  one  computer  code. 

A COM-GEOM  description  of  the  Air  Force  Man12  is  being  considered 
at  the  Ballistic  Research  Laboratory  (BRL).  If  completed,  this 
description  would  replace  the  currently  used  Computer  Man.  It  is  ex- 
pected that  this  new  Phantom,  unnamed  at  present,  will  provide  a signi- 
ficantly more  accurate  description  of  a representative  soldier. 

Much  work  is  still  needed  to  achieve  a working,  mortality- 
calculational  methodology.  Some  aspects  of  the  present  study  will 
probably  need  to  be  changed  or  refined.  The  author  I.opes  that  this 
study  will  contribute  to  the  achievement  of  this  goal. 


10. 

W . Saoao  and  L.  Sturdivan,  Unpublished  Data. 

H'W.  Kokinakis  and  W.  Cooper , Unpublished  Shook  Trauma  Assessments. 
12 

'Kase,  P. G. j " Computerized  Anatomical  Model  Man,"  Martin  Marietta 
Corporation  Technical  Report  No.  AFWL-TR-69-161,  January  1970. 

12'C.A.  Stanley  and  M.  Broun , "A  Computer  Man  Anatomical  Model," 
Ballistic  Research  laboratory  Report  ARBRL  TR  No.  02060 , May  1978. 
(AD  #A056564) 
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